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Abstract 

The nonlinear interactions between flexural and torsional modes of a microcantilever are ex- 
perimentally studied. The coupling is demonstrated by measuring the frequency response of one 
mode, which is sensitive to the motion of another resonance mode. The flexural-flexural, torsional- 
torsional and flexural-torsional modes are coupled due to nonlinearities, which affect the dynamics 
at high vibration amplitudes and cause the resonance frequency of one mode to depend on the am- 
plitude of the other modes. We also investigate the nonlinear dynamics of torsional modes, which 
cause a frequency stiffening of the response. By simultaneously driving another torsional mode in 
the nonlinear regime, the nonlinear response is tuned from stiffening to weakening. By balancing 
the positive and negative cubic nonlinearities a linear response is obtained for the strongly driven 
system. The nonlinear modal interactions play an important role in the dynamics of multi-mode 
scanning probe microscopes. 
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I. INTRODUCTION 



The Atomic Force Microscope (AFM) [l| is a crucial instrument in studying nanoscale 
objects. Various operation schemes are employed, which include the use of different 
cantilever geometries, higher modes or the torsional mode for imaging 2|-|5|. The nonlinear 
tip-sample interactions determine the dynamics in tapping-mode AFM and have been stud- 
ied in detail js, 7|. Besides this extrinsic nonlinearity, the intrinsic mechanical nonlinearities 
determine the dynamics of ultra-flexible microcantilevers at high amplitudes, as shown in a 
recent study Q]. These nonlinearities result in a amplitude-dependent resonance frequency 
and couple the vibration modes. In clamped-clamped beams, the nonlinear coupling is 
provided by the displacement-induced tension jo, |lo|. For cantilever beams it was shown 



that the coupling between the modes can be used to modify the resonance linewidth 111. In 

nn 

a multi-mode AFM |12l . Il3j , these modal interactions are of importance, since the resonance 
frequency of one mode depends on the amplitude of the other modes. 

In this work, we experimentally demonstrate the intrinsic mechanical coupling 
between the fiexural and torsional modes of a microcantilever. The resonance frequency of 
one mode depends on the amplitude of the other modes. The fiexural modes are coupled via 
the geometric and inertial nonlinearities. The torsional modes exhibit frequency stiffening 
at high amplitudes, which originates from torsion warping IJ]. Interestingly, the nonlin- 
earity constant of one torsional mode changes sign when another torsional mode is driven 
at high amplitudes. Finally, the coupling between the torsional and fiexural modes is studied. 



II. EXPERIMENT 

Microcantilevers are fabricated by photolithographic patterning of a thin low-pressure 
chemical vapor deposited silicon nitride (SiN) film. Subsequent reactive ion etching trans- 
fers the pattern to the SiN layer, and the cantilevers are released using a wet potassium 
hydroxide etch, resulting in a undercut-free cantilever. The dimensions are length x width 
X height (Lxwxh) = 42 x 8 x 0.07 fim^. These floppy cantilevers allow high amplitudes and 
thus facilitate the study of nonlinearities. The cantilever is mounted onto a piezo actuator, 
which is used to excite the cantilever. The cantilevers are placed in vacuum (pressure 
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< 10~^ mbar) to eliminate air-damping and to enable large vibration amplitudes, where 
nonlinear terms in the equation of motion dominate the dynamics. The cantilever motion is 
detected using a home-made optical deflection setup which resembles the detection scheme 
frequently used in scanning probe microscopes. The fiexural and torsional vibration modes 
are detected with a sensitivity of ± 1 pm/A/Hz 15[. A schematic of the measurement setup 
is shown in Fig. 1(a). The cantilever displacement signal is measured using either a net- 
work (NA) or spectrum analyzer (SA). To drive a second mode, a separate RF source is used. 




f,/kHz) f,,(kHz) 

FIG. 1: Measurement setup, (a) Optical deflection setup showing the laser beam, which reflects 
from the cantilever surface. The spot of the reflected laser beam is modulated in time by a 
frequency corresponding to the cantilever motion. The cantilever is mounted onto a piezo actuator 
in vacuum. Network (NA) and spectrum analysis (SA) is performed on the signal from the two- 
segment photodiode. (b) Frequency responses of the flrst and second flexural (top panels) and 
torsional (bottom panels) modes. Inset are the calculated mode shapes from Euler-Bernoulli beam 
theory. 



First, the flexural vibrations are characterized by measuring the cantilever frequency 
response at different resonance modes. The first flexural mode shown in Fig. 1(b) occurs at 
54.8 kHz with a Q-factor of 3000. The resonance frequency of the second mode is 347 kHz 
{Q = 3900), which is 6.33 times higher than the first resonance mode, in agreement with 
the calculated ratio fR,F2/fR,Fi = 6.27, following from Euler-Bernoulli beam theory. Not 
shown is the third flexural mode at 974.9 kHz, with /j?,F3//ij,Fi = 17.8, near the expected 
ratio of 17.6. This indicates that in the linear regime the cantilever beam is described by 
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the Euler-Bernoulli beam theory. Throughout the manuscript, the subscripts Fi and Ti 
indicate the frequency span around the i^^ flexural (F) or torsional (T) resonance mode. 
The subscript R refers to the resonance frequency of that particular mode. 

The torsional modes are characterized by rotating the cantilever over 90 degrees 
in the setup; the two-segment photodiode is then sensitive to vibrations corresponding 



to torsional resonance modes 



151]. The frequency response of the first two torsional 



modes is shown in Fig. 1(b). From theory, the ratio between the lowest two reso- 
nance frequencies of the torsional modes is 3, which is close to the measured ratio of 
/ii,T2 //ij,Ti = 1638kHz/535.4kHz = 3.06. The Q-factors of the first and second torsional 
mode are 4300 and 3200 respectively. 

At high drive amplitudes, the flexural and torsional modes become nonlinear. The 
nonlinearity of the flexural modes in a cantilever beam was theoretically studied by Crespo 
da Silva in 1978 ly, [iTj]. To include the torsional nonlinearity, the equations of motion 



are extended (Appendix A). For the flexural and torsional modes, the nonlinearity causes 
a Duffing-like frequnecy stiffening when the mode is strongly driven 18|, ll9| leading to 
a bistable vibration amplitude. This bifurcation is observed in all modes studied in this 
paper. These nonlinearities are responsible for the coupling between the flexural-flexural, 
torsional-torsional and fiexural-torsional modes. 



III. MODAL INTERACTIONS IN A MICROCANTILEVER 

We now experimentally demonstrate the coupling between the modes of a microcantilever. 
We use a two-frequency drive signal to excite two resonance modes of the cantilever simul- 
taneously while we measure the motion of one mode. First, we focus on the interactions 
between the flexural modes. Then we turn our attention to the torsional modes, starting 
with the amplitude-dependent resonance frequency of the torsional vibrations, followed by 
the demonstration of the coupling between the lowest two torsional modes. Finally, the 
interactions between flexural and torsional modes are discussed. 
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A. Flexural-flexural mode interaction 

To investigate the interactions between the two lowest flexural modes, the thermal 
motion of the first mode is measured with a spectrum analyzer, while the RF source 
strongly drives the second mode. The thermal noise spectra of the first mode as a function 
of the drive frequency of the second mode are shown in Fig. [2](a). The color scale represents 
the power spectral density of the displacement around the resonance frequency of the first 
mode. A shift of the resonance peak of the first mode is observed as the drive signal at 
/f2 approaches the nonlinear resonance of the second mode. The resonance frequency 
of the first mode for each drive frequency of the second mode is obtained by fitting 
the damped-driven harmonic oscillator (DDHO) response. In Fig. [2^b), this resonance 
frequency of the first mode is plotted versus the drive frequency of the second mode. The 
nonlinear response of the second mode is reflected in the resonance frequency of the first 
mode where the resonance frequency first increases and then jumps down after the second 
mode has reached its maximum amplitude, indicated by the arrow. At the maximum 
amplitude of the second mode, the resonance frequency of the first mode is shifted by 
several times its linewidth. This experiment shows that the coupling between the flexural 
modes can introduce significant resonance frequency shifts when multiple modes are excited 
simultaneously. Moreover, by measuring the shift in resonance frequency of the first mode 
the motion of the second mode can be detected. For comparison, the nonlinear response of 
the direct-driven second mode is shown in the inset of Fig. |2]^b). 



B. Torsional-torsional mode interaction 

Before turning to the interactions between the torsional vibration modes, we first 
measure the frequency response of a single torsional mode as a function of the drive 
strength. Although torsional modes are extensively used in AFMs [2I, Isil, their nonlinear 
behavior has not been investigated in detail. To investigate the nonlinearity, we strongly 
drive the torsional mode. In contrast to the fiexural-fiexural modes, where the nonlinearity 
arises from geometric and inertial effects, in torsional modes, the nonlinearity originates 
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FIG. 2: Flexural-flexural mode interactions, (a) Frequency spectra of the thermal motion of 
the first flexural mode (/fi)i when the second mode is driven through its resonance frequency. 
Colorscale represents the power spectral density of the displacement noise of mode 1. As the peak 
width remains constant, there is no significant change in the Q-factor. The motion of the second 
mode tunes the resonance frequency of the first mode, (b) The resonance frequency of the first 
mode /r,fi versus the drive frequency of the second flexural mode. The nonlinear response of 
the second mode is reflected in the fitted resonance frequency of the first mode. Inset: the direct 
measurement of the nonlinear frequency response of the second mode 2^ 

from torsion warping and inertial moments [l^. In Appendix A we discuss the equations 
of motion including the nonlinearities involved. The amplitude of the first torsional 
mode with varying drive power is shown in Fig. |3]^a), with selected frequency responses 
of the first torsional mode plotted in Fig. El^b). At low driving power, the resonance 
line shape is a DDHO response, and the cantilever is oscillating in the linear regime. 
When the power is increased, the frequency response leans towards higher frequencies 
and the amplitude bifurcates. Close to this critical amplitude (0 dbm) the slope of the 
frequency response approaches infinity, which may be used to enhance the sensitivity in tor- 
sional mode AFM. A frequency stiffening is observed for the first and second torsional mode. 



Similar to the flexural-flexural modal interactions discussed in the previous section, 
the coupling between the first and second torsional modes is studied: we measure the 
thermal noise of the first mode while the drive power of the second torsional mode is varied. 
The resonance frequencies, obtained from DDHO fits to the thermal noise spectra of the 
first mode, are shown in Fig. Hl^a). The resonance frequency increases with 500 Hz, while 
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FIG. 3: Nonlinear torsional mode, (a) Frequency responses of the first torsional mode, when the 
drive amplitude is increased. Beyond a power of 5 dBm, the response is bistable. Color scale 
indicates the amplitude normalized to the drive voltage, (b) Resonator amplitude traces at 5 
selected drive powers. The nonlinear frequency response is visible at high drive powers. The 
frequency is swept from low to high. 

increasing the driving strength of the second mode to 10 dBm. We now perform a similar 
experiment as the one as shown in Fig. [2]^b). Thus, the first torsional mode is used to 
detect the nonlinear vibrations of the second torsional mode. Fig. ID^b)) shows the nonlinear 
response, resembling the behavior of the first mode shown in Fig. [3t^b). 

Interesting behavior is observed when both torsional modes are driven in the nonlinear 
regime. In contrast to measurements in the previous section, the first mode is now also 
driven in the nonlinear regime. Fig. |5]^a) shows the nonlinear frequency response of the first 
torsional mode, while stepping the drive frequency of the second torsional mode through 
its resonance. Fig. Mjo) shows individual traces, which reveal interesting behavior; in the 
lowest panel (i), there is no influence of the second mode and frequency stiffening of the first 
mode is observed cf. Fig. Mjo)- When the amplitude of the second mode starts to increase 
as it approaches its resonance, the response of mode 1 becomes more linear (panel ii). 
Here, the frequency stiffening and weakening nonlinearities are balanced yielding a linear 
response. At high amplitude of the second mode, frequency weakening of mode 1 (panel iii) 
is observed. When the amplitude of the second mode drops, frequency stiffening is restored 
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FIG. 4: Torsional-torsional mode interactions, (a) The resonance frequency of the first torsional 
mode, when the drive power at the second torsional mode on resonance is varied. The resonance 
frequency of mode 1, /r,tii increases with the drive power of the second mode, pT2- (b) The 
nonlinear response of the second mode is measured by using the first mode as a detector. 

(panel iv). This measurement not only demonstrates the coupling between the torsional 
modes, but also that the sign of the nonlinearity constant of a torsional mode depends on 
the amplitude of the motion of the other modes. By simultaneous driving another mode, 
the torsional frequency response can be tuned from a stiffening to a weakening characteristic. 



C. Flexural-torsional mode interaction 

The coupling between the first flexural and first torsional mode is now studied experi- 
mentally. This coupling is theoretically described in Eq. IA5I Fig. [6]^a) shows the resonance 
frequency of the first torsional mode as a function of drive power of the first flexural 
mode. The resonance frequency increases with 100 Hz when the power in the flexural 
mode is increased to 10 dBm. Detection of the nonlinear flexural mode by measuring the 
resonance frequency of the torsional mode is shown in Fig. Mjo). The nonlinear interactions 
when both modes are driven in the nonlinear regime are shown in Fig. [6](c) and (d). The 
interaction is clearly visible in the frequency-frequency plots, where one frequency is swept 
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FIG. 5: Tuning the torsional nonlinearity via modal interactions, (a) Frequency responses of the 
nonlinear first torsional mode, while the frequency of the second mode is swept trough its nonlinear 
resonance. The frequency stiffening of the torsional mode changes into weakening when the second 
mode oscillates with high amplitudes. When the amplitude of the second mode jumps down, again 
frequency stiffening is observed, (b) Traces from (a) taken at the indicated frequencies. In panel (ii) 
the response is close to a linear one, due to balancing of the stiffening and weakening nonlinearities. 

and the frequency of the RF source is stepped across the nonlinear resonances of both modes. 



IV. DISCUSSION AND CONCLUSION 

In summary, we demonstrated the coupling between the flexural and torsional vibration 
modes in a microcantilever. This coupling is due to nonlinearities, which also give rise to 
a amplitude-dependent resonance frequency. The interactions between the different flexural 
modes, between different torsional modes and between the flexural and torsional modes 
are demonstrated in detailed experiments. We also demonstrate the nonlinear frequency 
stiffening of torsional modes driven at high amplitudes. 

Several applications are proposed for the modal interactions. A specific resonance 
mode can be shifted to a higher frequency by simultaneously driving another mode. For 
strongly driven modes, the cubic spring constant (nonlinearity) can be modified from positive 
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FIG. 6: Flexural-torsional interaction, (a) Resonance frequency shift of the first torsional mode, 
when the drive power of the first flexural mode is increased. The resonance frequencies are obtained 
from thermal noise spectra, (b) The nonlinear resonance response of the first flexural mode reflected 
in the resonance frequency of the torsional mode (from thermal noise spectra) . (c) and (d) Nonlinear 
dynamics when the flexural as the torsional mode are both excited in the nonlinear regime. 

to negative, tuning the response from stiffening to weakening. By balancing two excitation 
strengths, a nonlinear response can be tuned to a linear one. By modal interactions, one 
mode can be used to detect the motion of another mode of the same cantilever. Besides 
these applications, the modal interactions have consequences for multi-mode schemes, such 



as the scanning probe microscopy and mass sensors based on microcantilevers 13|, |22|-|24|. 
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Appendix A: Theory of modal interactions 



In this appendix the nonhnear equations of motion of the modes in a cantilever are de- 
scribed. We start with the general equations of motion, which include the coupling between 
the torsional and flexural modes. Then, we consider flexural modes along one axis. We 
conclude with the equation of motion of two coupled flexural modes, relevant to the experi- 
ment in Section 3.1. We start with the equations derived by Crespo da Silva [igI, Q, who 
described the motion in two flexural directions v and w and the torsional angle 9: 



+ {D.^ - D(;)WsVssWss + 
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ds 



d 



D,-D^^ — (9wss-9^'vss 
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d_ 
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d 



St 



ds . 



D^{9s + VssWs) +{Dr,-Di;) {Vss - Wss)9 - VssW, 



St 



Qe 



(A3) 



Here, the subscripts s and t denote differentiating to time and position respectively. ^, rj and 
C represent the principle axes of the beam's cross section. 7 and Qe represents the damping, 
D,^^(^ are the flexural stiffnesses of the beam and is the torsional stiffness. The moments 
of inertia are given by The driving force is q with driving frequency VL. 

Now, we consider only vibrations in one direction, so the Eq.A.2 and all terms with w in 
Eqs.A.l and A. 3 are disregarded. For the torsional mode, the nonlinear effect of torsion 
warping is taken into account [l4|, where we assume that the ends of the beam are axially 
immovable. The coupled differential equations (in non-dimensional form, notation from 
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Ref. [25|) are now given by: 

Vtt + + (3yVssSS = 

/S s' 
ds' m ds" —{vl))+q cos(fit), (A4) 



d_ 

ds 



(1 - (3y)e^vss 



d_ 

dt 




-(l-'3„)i49+(j,-Jc)«i9 = 0. 



dsstt + 



EC 



(A5) 



/3y and /3e are the ratios between two stiffnesses {(3y = Dc^/D^^ and Pe = D^/Dr^) and A is 
the cross sectional area. The torsion nonhnearity is written as (3z = ^Eln- The torsion 
constant It, the warping constant Cs and the time-dependent tensile axial load are given 
by: 



It 

Cs 



y dz 



dy 



dS, 



2 / 



y \2 



dx. 



(A6) 



Here, 5* is the solid angle and cj)^ is the primary warping function. A more detailed descrip- 



tion of the nonhnearity in the torsional mode is found in Ref. 14 1. 

To demonstrate the origin of the nonlinear interactions observed in the main text, we 
now simplify the coupled equations Eq.A.4 and A. 5 by applying the Galerkin procedure. The 
solutions are then written as a linear combination of the linear mode shapes of the cantilever 
with coefficients, which correspond to the time-dependent vibration, v = J2i {s)vi{t) and 
6 = '^iE^6i{t), where i represents the mode number. The mode shapes of the fiexural and 
torsional modes of the cantilever will be discussed. Introducing the operator C, the linear 
part of Eq. IA5I is written as: 



2 TPV 



(A7) 



The resonance frequency is denoted as u. The eigenfunctions can be calculated together with 



the boundary conditions for a single-clamped cantilever F'"{0) = 
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as: 



K - 

k\ sinh(A;i) 

The values of ki^2 are given by: 

kio = 



F = [cosh(A;is) — cos(A;2s) — K{smh.{kis) 
kf cosh(A;i) + k^ cos(A;2) 



ki/k2 sin(A;2s)], 



kik2 sin(A;2 



(A8) 



^ ^ 2/3, ^ ^ ( 2A 



+ 



B 



(A9) 



The values of ki and ^2 depend on the mode number i and can be calculated via the 
generating function: 



ki + k2+ 2k^k2 cosh(A;i) cos(fc2) + fciA;2(^2 — ^i) sinh(A;i)sm(A;2) = 



(AlO) 



The dimensional resonance frequency of the flexural mode is given by ujB,i = 
K,i(h/L'^)^y D(^/ where Kj is 1.875, 4.695 and 7.855 for i = 1, 2 and 3. The beam shape of 
the first two flexural modes are shown in the inset of Fig. [T]of the main text. 

The torsional mode shapes can be calculated by considering the operator with 
eigenvalues ut 



M[G] 



— — UJrpLr, 



3i dt 



(All) 



and the corresponding boundary conditions of -F^(O) = F^{1) = 0. Inserting the boundary 
conditions in Eq. lAllt gives the equation for the torsional mode shapes: G = sin[(22 — 
l)7r/2s]. The resonance frequency of the torsion mode is given by cut = (2i — l)(7r/2)A//3g/j^. 

The Galerkin procedure is applied to Eq.A.4 and A. 5: i.e. the solutions are written as 
a linear combination of the eigenmodes. We assume that the flexural mode is only excited 
around the resonance frequency, accumulating in the equations: 



4 + i^vi + ujW 



EEE 

j k I 



a2,ijkiv'"v'' + a^^ijkiiv^v^A + f cos(fit). 



(A12) 



j k I 

The above equations show that the nonlinearity is the origin of the modal interactions. 
Note that for j = A; = /, the nonlinear equation describing one mode of the cantilever is 
found. A quadratic couphng is present between two different vibrational modes (for example 
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TABLE I: The values of the integrals in the coefficients a for the interactions between the first and 
second flexural and torsional modes. 

k — I) also follows directly from the cubic nonlinearities. This quadratic coupling is clearly 
observed in the experiments. In Eq. A. 12, the terms linear in 9 are assumed to only modify 
the resonance frequency out- The coupling (Galerkin) coefficients a are given by the following 



equations: 
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Considering the interactions only between the lowest two modes of the torsional and the 
flexural mode, the values of the integrals in the coefficients a are given in Table A.l. 
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To give an example, we work out Eq. A. 12 for the fundamental and second flexural mode 
of a cantilever. This case is studied in the experiment and described in Section 3.1. We 
denote the amplitudes with — a and — b) and the coupled equations are given by: 

att + 7i«t + ^f(^ = Oi2,una^ + Q;2,i222&^ + (q;2,iii2 + a2,ii2i + Oi2,i2n)a% 

+ (a2,i2i2 + "2,1221 + a2,ii22)&^a + a3,iiii2a(a^)tt + 03,122226(6^)^ + (03,1112 + 03,1121)0(06)^ 
+ 03, 1212)6(a6)tt + 03,i2ii6(a2)tt + 03,11220(6^)^ + /i cos(Ot), (A14) 
btt + 726* + ujpb = a2,2iiia + 02,22226 + (02,2112 + 02,2121 + 02,2211)0 6 

+ (02,2212 + 02,2221 + 02,2122)6^0 + 03,2iii2a(a^)tt + 03,2222 26(6^)tf + (03,2112 + 03,2121)0(06)^ 

+ (Qi3,2221 + Qi3,2212)6(o6)tt + 03,22116(0^)^ + O3,2122o(6^)tt + COs(Oi). (Al5) 
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